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Abstract 

, We study a class of inhomogeneous and anisotropic G2 string cosmological models. In the 

^ case of separable G2 models we show that the governing equations reduce to a system of ordinary 

■ differential equations. We focus on a class of separable G2 M-theory cosmological models, and study 

O '■ 

Q"^ . their qualitative behaviour (a class of models with time-reversed dynamics is also possible). We find 

o , 

CN that generically these inhomogeneous M-theory cosmologies evolve from a spatially inhomogeneous 

o ■ 

and negatively curved model with a non-trivial form field towards spatially flat and spatially 
homogeneous dilaton-mo dull- vacuum solutions with trivial form-fields. The late time behaviour is 
r~| , the same as that of spatially homogeneous models previously studied. However, the inhomogeneities 



X 
J3 



are not dynamically insignificant at early times in these models. 
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I. STRING COSMOLOGY 



Non-perturbative M-theory encompasses and unifies all five anomaly-free, perturbative 
superstring theories 0] and corresponds to eleven-dimensional supergravity in the low- 
energy limit 10]. In particular, the compactification of M-theory on a circle, S^, leads to 
the type IIA superstring. A study of the qualitative cosmological effects that can arise in 
M-theory is therefore of considerable importance. To lowest order (in the inverse string 
tension), the tree- level effective action for massless fields contains a dilaton, a form field 
(which in 4-dimensions is dynamically dual to pseudoscalar axion field) and (a) stringy cos- 
mological constant (s). Even in this approximation the one-loop string equations of motion 
for inhomogeneous backgrounds are very difficult to solve, and it is a useful first step to 
consider models in which the homogeneity is broken only in one spatial direction. Metrics 
that admit two commuting (orthogonally-transitive) space-like Killing vectors are referred 
to as G2 spacetimes. 

String models admitting an abelian group, G2, of isometries have a number of important 
physical applications. The spatially homogeneous Bianchi types l-Vllh and locally rotation- 
ally symmetric (LRS) types VIII and IX admit a G2 group of isometries and so the G2 
cosmologies can be considered as inhomogeneous generalizations of these Bianchi models. 
Non-linear inhomogeneities in the dilaton and axion fields can be investigated and, in prin- 
ciple, this allows density perturbations in string-inspired inflationary models such as the 
pre-big bang scenario to be studied 0, |^. Given the potential importance of this scenario 
it is important to study its generality with respect to inhomogeneities as well with respect 
to anisotropics. The general effects of small inhomogeneities and anisotropics have been 
studied by Veneziano |Q. 

In general relativity (GR) the generic singularity is neither spatially homogeneous or 
isotropic. Hence it is of interest to study more general models. In particular, it has been 
conjectured that G2 metrics represent a first approximation to the general solution of Ein- 
stein gravity in the vicinity of a curvature singularity 0, ^ |§, ^ . The high curvature regime 
is precisely the regime where stringy deviations from GR are expected to be significant. 
The G2 models studied here may therefore provide insight into the generic behaviour of 
cosmologies at very early times. 

A number of exact inhomogeneous and anisotropic G2 string cosmologies have been found. 
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Barrow and Kunze studied an inhomogeneous generalization of the Bianchi type I string 
cosmology |]T0[ and Feinstein, Lazkoz and Vazquez-Mozo derived a closed, inhomogeneous 



model by applying duality transformations on the LRS Bianchi type IX cosmology [ITT 



Clancy et al. have found inhomogeneous generalizations of the Bianchi type VI/i universe 
and have studied their asymptotic behaviour . 

In general, the field equations reduce to a system of coupled, partial differential equations 
in two variables when spatial homogeneity is broken along a single direction. Unfortunately, 
these equations are still very complicated. However, solutions can be found due to the 
non-compact global symmetries of the string effective action. When the metric admits two 
commuting space-like Killing vectors, there exists an infinite-dimensional symmetry on the 
space of solutions that may be identified infinitesimally with the 0(2,2) current algebra 



p!3| , [T^ , [T5[| . This symmetry reduces to the Geroch group, corresponding to the SL(2,i?) 



current algebra, when the dilaton and two-form potential are trivial |T6|, and includes the 
global SL(2, R) S-duality of the action. 

New inhomogeneous G2 string cosmologies containing a non-trivial two-form potential 
may be generated by an application of both the S- and T-duality symmetries from sim- 



pler (dilaton- vacuum) seed solutions. Lidsey et al. [0 discuss the non-compact, global 
symmetries of the string effective action in a variety of settings and review various methods 
for solving the Einstein-scalar field equations utilizing generating techniques (from solu- 
tions with a minimally coupled, massless scalar field from a vacuum, G2 cosmology). In 
particular, Feinstein, Lazkoz and Vazquez-Mozo [ill] present an algorithm which permits 
the construction of inhomogeneous string solutions by employing a Buscher transformation, 
inverse scattering techniques 0, |^, followed by the generating technique of Wainwright, Ince 
and Marshman [Q. Feinstein et al. employ this algorithm to generate a closed, inhomo- 
geneous string cosmology with topology from a LRS Bianchi type IX solution O, |19 



However, this algorithm involves a number of non-trivial operations, and an alternative and 
more straightforward approach is to apply an 0(2,2) transformation directly to the seed 



cosmology [17 



In this paper we shall consider a class of separable string cosmological models whose 
governing equations reduce to ordinary differential equations (ODE) which can be studied 
by qualitative methods. In particular, we shall focus on a class of M-theory cosmological 
models. 
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II. STRING ACTION 



We consider the general string action in the form 

1 



S 



d Xy/—g 



R + {V^y - 6(V/5)2 - -e^^iVay 



2A 



M 



(2.1) 



in terms of the pseudo-scalar axion field, a, the 4-D dilaton field and the modulus field 
/?, where A and Am represent cosmological constant terms and may be interpreted as 
a 0-form field strength. This is a phenomenological action representing the bosonic sector 
of the effective supergravity action for the low-energy limit of M-theory and encompasses 
other string theories EO] . We are particularly interested in the class of four-dimensional cos- 
mologies derived from the type IIA string and M-theory effective actions and which include 
a non-trivial Ramond-Ramond (RR) sector . In these models a specific compactification 
from eleven to four dimensions was considered, where the topology of the internal dimensions 



was assumed to be a product space consisting of a circle and an isotropic six-torus ; this 



is dynamically equivalent to compactifications on a Calabi-Yau three-fold |23|. The FRW 



models in this class of cosmologies was studied in pi 
Defining 



/:M = -3e"*(V/5f --e*(Va) 



Ae"* - ig^e-^^ - ^Am; 



2 

the Euler-Lagrange equations then lead to the field equations (FE) 

1 



Tal3 = 9al3^M 



.dC 



M 



dg 



a/3 ' 



G 



+AMe'^ 



20$ 



□a 



(2.2a) 

(2.2b) 

(2.2c) 
(2.2d) 



In the above, greek indices take on values 0, 1, 2, 3, and units are chosen so that 16ttG = 1. 
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III. G2 COSMOLOGIES 



Let us examine (|2.2|) within the context of G2 cosmological models described by the hne 
element 

ds^ = e^^ {-de + dz^) + [e'^dx^ + e-^dy^) , (3.1) 

where the metric functions {F, G, g} and the string functions {$,/?, a} are all functions 

Qq Qq 

of t and z only. For any q(t,z), we define q = — , q = — , and A^g = 9 — q"- The 

at 

local behaviour of these models is determined by the gradient 5^ = S^G, and cosmological 
solutions arise if -B^ is globally time-like. 
Also, the Ricci scalar is given by 

1 



R = 
2 



Using these expressions, and defining the modified dilaton field, 

^ = ^ - F -G, (3.3) 
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the field equations become 
AV = ] 



+3 (^^2 _ j _ 1 g2^+2F+2G (^2 _ ^/2) _ ^^2F ^ 



AV 



+ F + 26- ) a - ((^' + F' + 2G") a' 



A^G 



A^g = l^ + F)q-{^' + F')q' 



A^F 



-3 - + ^e2^+2F+2G (^^2 _ ^/2) _ ^g2F 



2 



^2 _j_ ^'2^, g2.p+2F+2G 



+ = ^GG' + ^qq' + Q/3p' + F'(ip + F)+F{ip' + F') 

^Ig2^+2F+2G^^/^ 



These equations reduce to those in |20, 21 1 in the appropriate limits 
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IV. SEPARABLE G2 STRING COSMOLOGIES 



A. General Case 

Let us assume separability of the metric functions of the form 

F(t,z) ^ F{t) + f{z), 
G{t,z) = Git)+g{z), 
q{t,z) = q{t) +v{z), 

and appropriate separabiUty conditions on the matter fields $(/!:, z), z), a{t, z). Then the 
Ricci Scalar is given by 

R = ^e-^^-^f[AG + AF + SG^ + f - (Ag" + Af" + ?,g'^ + 

If 

4/ + 4/" + 3(7" + z/'' = C, (4.1) 

where C is a constant, then we obtain a condition which constrains the spatial dependence 
of the metric. The Ricci Scalar is then given by 

R = Je-2^-'^[4G + 4F + 3G^ + f - C]. 

Putting this expression for the Ricci Scalar into the action ( |2.1| ), the spatial dependence 
of the geometrical terms can be eliminated (by integration over the spatial coordinates in 
the action). After applying any further separability conditions (on the matter fields), the 
resulting FE will be a system of ODEs. Note that the effect of the spatial dependence is to 
add a further contribution (C) to the cosmological constant A in the action. 
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B. Specific Example: Linear Dependence in z 



In an attempt to remove the 2;-dependence, let us assume separability of the form 



F{t,z) 


= m + ^cz 


G{t, z) 


= Git), 


q{t, z) 


= q{t) + az, 


Ht,z) 


= + mz, 


f3{t,z) 


= (3{t) + nz, 


a{t, z) 


= (j{t) + Iz, 



where a,c,l,m,n are constants (a^ is equivalent to the constant C in equation ( [4.1|) ), and 
therefore 

(p{t, z) = - F{t) - G{t) + (^m-^cj z = ifit) +(^- z. 
With the above assumptions, the metric becomes an extension of the inhomogeneous scalar- 



field G2 solutions found by Feinstein and Ibanez ||22[ to M-Theoretical models. In addition, 
for particular values of the parameters, the metric reduces to spatially homogeneous Bianchi 
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/, III and VIq models. Hence, 

_^^2F+cz _ 1 + 2m^ + 12n^) , (4.2a) 

/3 = (^0 + ^) + ig2g[<^+3F+G-6/3+{c+m-6n)^] _ ^^ ^b) 

a = -(^(f + F + 2G^a + ml, (4.2c) 
q = (^(^ + q- ma, (4. 2d) 

ir = 1 (^0 + _ 1^2 _ 1 -2 _ 3^2 ^ lg2^+2F+2G+2„^. (^2 _ 

+^(a'-2m' + 12n2), (4.2e) 

G = (0 + f) G + - /2) e2^+2F+2G+2„^. 

+ = 2(F + + m) (^F + lc^ +1^2 

+ ^ (g + a)' + 6 (/? + n)' + i (a + /)' e2^+2F+2G+2™.^ 

(4.2g) 

Note that the constraint equation can be rewritten as: 

= (^F + -^^ + c(^F + ^ + rr?j +^G^ + ^{q + af 
+6 + n) ' + ^ (d + If e2^+2F+2G+2m. 

^^g2F+c2; _j_ lg2g[v>+3F+G-6/3+(c+m-6n)2] _^ ^^g</p+3F+G+(c+m)z ^ (^^^3^ 

In order for the FE to be independent of 2, it is necessary that m = and that either 
c = or c = 6n. Furthermore if c = then we have that either n = or Q = 0. In the 
c = 6n case, we have that A = Am = 0. It is the latter case that is of interest to us here. 
From here forth we shall assume that c = 6n and that A = Am = 0. This particular subcase, 
which is of relevence in M-theory cosmology, is of special physical interest. (The resulting 



FE in the remaining cases are displayed in ||24|| ) 
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V. INHOMOGENEOUS M-THEORY COSMOLOGICAL MODELS 



Substituting c = 6n, m = 0, A = 0, Aj\/ = into ( p.4|) , ([4.3|) (and taking the linear com- 
bination [( 4.2g )-( 4.2e|) -( ^?2a )]) we obtain the following system of ODE with two constraints: 



(p = ^(q'^ + G^-a^ + 2(0 + + 12/3^ - 12n^ 

a = -a{ip + F + 2G) (5.1c) 

q = q{ip + F) (5.1d) 
F = 1 (3a^ + f + G' + 36n' + 12/3^ + e'^^+'^^+^G {f + 3&') 

-2(0 - F)^ + SF^) (5.1e) 



= 202 - - g2 - 12/?2 - 2F^ - 12n^ - 

_g2v,+2F+2G(;2 ^ ^2) _ g2g-6/3+^+3F+G (g^^g) 

= ag + /e2^+2^+2G^ + 6n(0 + F + 2/3) (5.1h) 



From equation (|5.1g|) we are able solve for and make a global substitution for the quantity 



Q2g 6/3+1/J+3F+G Malting this substitution we have the following system of ODEs: 

y3 = i (^g2 ^ ^2 _ ^2 ^ 2(0 + Ff + V2i? - I2n^ + e2^+2^+2G(/2 _ ^2^ j (5 2a) 

/3 = i (^202 - G2 _ ^2 _ ^2^2 _ 2^2 _ ^2^2 _ ^2 _ g2^+2F+2G(^2 ^ ^2^ ^ 

+/?(0 + F) (5.2b) 

a = _a(0 + F + 2G) (5.2c) 

q = g(0 + F) (5.2d) 
F = i (^3a2 + g2 + (^2 ^ 35^2 ^ ^2/5^ + e2^+2^+2G(/2 + 3^2^ _ 2(^ _ pf + gF^j^^ ^e) 

G = G(0 + F)-e2^+2F+2G(^2_^2) 

_1 (^202 - G2 _ ^2 _ ^2^2 _ 2^2 _ ^2^2 _ ^2 _ g2^+2F+2G(^2 ^ ^2^ ^ (g ^f) 

= aq + le^'P+^^^^^a + Qn{^ + F + 2(5) (5.2g) 
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From the constraint (|5.1g|) we see that if = 0, then all of the other state variables must 
be simultaneously zero, which can only occur at an equilibrium point of the system. Hence 
if must be positive (or negative) throughout the physical phase space. Here we shall assume 
> (the case < can be obtained by a time reversal - see later). 

We define new variables of the form 



/3 = 76^ 



G 



1 G 



V2 



^1 



a2 + 12n2 1 



^2 = ^e^+^+^- 
V2 <f 



(5.3) 



and a new time variable 

dt _ 1 

dr (fi 

The variables are chosen so that the transformed dynamical system has a compactified 
phase space. This property comes from the fact that Q2g--6/3+</3+3F+G y g -^vj^iQ^^ implies that 
equation ( |5.1g| ) yields 



l> + G^ + f + + + ^ + ^ 



The dynamical system ( ^.2|) becomes 



dF 

dr 

dG 

dr 

dq 

dr 
da 

dr 

dp 

dr 

d-^i 

'd^ 

d^2 

dr 



F(2F-f) + i(g2 



G' + (3' + 3a' + 3^f + - (1 - FY 



G[l + F-f)- — [l-F^-G^-q 
q{l + F-f), 
-aiV2G + r), 



3a' - /3' - §^ + 



2 ' 



/3(l + F-r~) + ^(l 
^oil + F + V2G-r), 



G' 



(5.4) 

(5.5a) 

(5.5b) 

(5.5c) 
(5.5d) 

(5.5e) 

(5.5f) 

(5.5g) 



and 



1 + FY + G'' + f + P'^ + a 



:i2 



where the constraint equation becomes 
0= \l . '.^ . ^1 



a2 + 12n2 



V2aq + 6n{l + F + ^P) 
v6 . 



+ 2*2^. 



(5.6) 
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There exists a first integral in the physical phase space (g 7^ 0, 7^ 0) for this system. 
The function 

M = 

is constant , i.e., M' = 0. This implies a first integral for the original system of ordinary 
differential equations ( p.2| ) 

where C is a constant. 



1. Invariant Sets, Monotonic Functions 

We first recall that the phase space for this dynamical system is the interior and boundary 
of the compact set given by 

1> + + f + + 13^ + ^ + (5.7) 

Various hyperplanes divide the phase space into a number of different regions, they are 
g = 0, a = 0, \&i = and \E'2 = hyperplanes. We note that = divides the phase space 
into two distinct regions < and > 0. The dynamics in the invariant set < is 
the time reversal of the dynamics in the invariant set \l'i > (see 
Consider the function 

Ml - ^ 



and its derivative 

dMi 



-2Mi(l + F) 



dr 

We easily see that this function is monotonically decreasing in the invariant set g 7^ 0, 5" 7^ 
0, \l/2 7^ 0. Therefore, we can conclude that there are no closed or periodic orbits in the 
seven dimensional phase space, except possibly on the lower dimensional boundaries of this 
seven dimensional invariant set. 

We restrict ourselves now to the invariant set a = 0. Consider the function 

^2 

and its derivative 

dM2 

= -2M2 1 + F) 

dr 
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We easily see that this function is monotonically decreasing in the invariant set g 7^ 0, 7^ 0. 
Therefore, we can conclude that there are no closed or periodic orbits in this six dimensional 
phase space, except possibly on the lower dimensional boundaries of this six dimensional 
invariant set. 

In the six dimensional invariant set \1'2 = 0, the function 
has the derivative 

— ^ = -2M3 1 + F 
ctr 

which is monotonically decreasing in the set \&2 = 0, g 7^ 0, \l'i 7^ 0. Therefore we conclude 
that there are no closed or periodic orbits in this six- dimensional invariant set. 
In the six dimensional invariant set q = 0, the function 

M4 = 

has the derivative 

^ = 2M,(1 + F) 

CLT 

which is monotonically increasing in the setg = 0,o"7^0,\I'27^0. Therefore we can conclude 
that there are no closed or periodic orbits in this six dimensional set. 

With the existence of these monotonic functions Mi, M2, M^, and M4, we can conclude 
that there are no closed or periodic orbits in the physical six dimensional phase space, [except 
possibly on lower dimensional (less than 5) invariant sets]. 

The zero- curvature spatially homogeneous and isotropic spacetimes are contained in the 
set q = ^1 = ^2 = G- V2F = union g = ^i = a = G- V2F = 0. 

The matter fields in satisfy various energy conditions. For example, the positivity 
of the kinetic energy of the pseudo-axion scalar field, a, demands that — > (i.e., 
(T^ — ^2 — 0)- However, we note that — \t'2 = is not an invariant set. 



2. Equilibrium Points and Exact Solutions 

There are two equilibrium points and one three-dimensional equilibrium set. 
• The three-dimensional equilibrium set is given by 

1^2 ^ ^ q^ + = l,a = 0,^1 = 0,^2 = 0} 
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Note that since both = and \l/2 = we necessarily have that a = n = I = 0. At this 
point the value of f = 1 + Fq. The exact solution is then 

(pit) = hi ~-^\n{ft + ho), 
r 

F{t) = Foihi-l\n{H + ho)) + Fi 
r 



Git) = V2Go{hi-^Hft + ho)) + Gi 
qit) = V2qo{hi - I ln(ft + ho)) + 



cr{t) = cTi 

Pit) = -^Po{hi-l\n{rt + ho))+Pi 

where Fq + Gq + +/5q = 1 and where Fi,Gi,qi,ai, Pi, hi and ho are integration constants. 

Since a = n = I = 0, this metric is spatially homogeneous (and flat). Since a = and 
Q = (which follows from the other conditions), this equilibrium set represents spatially 
fiat solutions where the form-fields (the axion field and the four-form field strength) are 
trivial and only the dilaton and moduli fields are dynamically important. These solutions 
are known as the 'dilaton-moduli-vacuum' solutions (and their analytical form is given in 

Recall that the dynamics of these models is restricted by the constraint given by equation 
( p.6| ). At these equilibrium points we are able to locally solve for the value of and 
substitute into the remaining equations. The eigenvalues in the six-dimensional constraint 
surface are 

0, 0, 0, V2G, -V2G - 1 - F, -1 + F + V2G - Vq(3 

The three zero eigenvalues correspond to the fact that this equilibrium set is three- 
dimensional. It is clear that a subset of this equilibrium set will act as saddles of varying 
degree of stabihty, while another subset will act as sinks. Consequently, a subset of these 
dilaton-moduli-vacuum solutions with trivial form fields are sinks in the physical phase 
space, even in the presence of (negative) spatial curvature, and are thus generic attracting 
solutions. 

• A second equilibrium point is 

{F = -\.G = 0,q = 0,P = -^,^ = ±^'^1 = = 0} 



14 



Note that since \l/2 = we have necessarily that / = 0. At this point the value of f = 0. 
[Note, since the dynamical system ( p.5| ) is invariant under the transformation (\E'i,\&2) — *■ 
(— — \^2) there exists a corresponding equilibrium point with a = — v^/3.] The exact 
solution is then 

(p{t) = hot + hi, 

= -^ihot + hi) + Fi 
G{t) = Gi 
(lit) = qi 

a{t) = ^e~imhot+h,)-(FM) ^ 
Pit) = -hhot + hi) + (3i 



where ho = fVa^ + I2n^, = ^{a^ + 12n'^)e'^^'~'^^'~^^ and Fi,Gi, qi, ai, Pi, hi are all 
constants. In this situation the variable q can be eliminated. The eigenvalues restricted to 
the constraint surface are 



1 1 1 



3, 3, g(l ± V 15 + 8v^z), i(l ± V 15 - 8v^z) 

This point represents a past attractor or a source. This corresponds to a spatially non- 
vacuum inhomogeneous model with a diagonal Einstein Tensor, having negative curvature. 
The line element corresponding to this solution (after a few coordinate redefinitions) is 

• The third equilibrium point is 

{F = -5.G . ^.5 = J = a = 0, = H^. = 0} 

Note that since \1'2 = we have necessarily that / = 0. At this point the value of f = 0. 
[Note, since the dynamical system ( ^.5|) is invariant under the transformation (\^i,\^2) 
(— ^1, —^2) there exists a corresponding equilibrium point with a ^1 = — 2-\/3/7.] The exact 
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solution is then 



(f{t) = hot + hi, 



Fit) = ~ihot + hi)+Fi 

G{t) = '^{hot + hi)+Gi 
q(t) = qi 



a{t) = ai 

Pit) = —{hot + hi) + pi 



where Hq = j^Va^ + l^n^, = + 12n^)e'^^'-^^'-^^ and Fi,Gi,qi,ai, l3i,hi are all 



constants. This solution is a curved inhomogeneous model with a trivial axion field. In 
this situation the variable q can be eliminated. The eigenvalues restricted to the constraint 
surface are 



This point is always a saddle. 
VI. DISCUSSION 

We have studied several classes of inhomogeneous string models whose governing equa- 
tions reduce to ODE. In particular, we have found that generically solutions of the class of 
separable G2 inhomogeneous M-theory cosmologies studied evolve from a spatially inhomo- 
geneous and negatively curved model with a non-trivial form fields towards (a subset) of 
spatially flat and spatially homogeneous dilaton-moduli-vacuum solutions where the form- 
fields (the axion field and the four-form field strength) are trivial and only the dilaton and 
moduli fields are dynamically important. This late time behaviour is the same as that of 
the spatially homogeneous models studied previously. However, in these models the inho- 
mogeneities are not dynamically insignificant at early times, and the models asymptote (in 
the past) toward a new class of inhomogeneous cosmological models. 

As noted earlier, the time-reversed dynamics of the (p > models we have considered 
thus far is equivalent to the dynamics of the case where (f < 0. This follows by redefining the 
time variable according to |^ = — 4 and appropriate definitions of the other state variables. 
The evolution equations will have an 'overall' change in sign, and hence the equilibrium 



2 2 4 




7' 7' 7' 
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points are identical in both cases, but the eigenvalues have opposite signs. Consequently, 
the dynamics of the < models is the time reversal of the > models and the time- 
reversed dynamics of the above class of models is deduced by interchanging the sources and 
sinks and reinterpreting expanding solutions in terms of contracting ones, and vice-versa. 

Although at late times (in the > models) the inhomogeneities decay, the inhomo- 
geneities are important at intermediate times and, in particular, at early times. Thus the 
qualitative features of the models are quite different to those of spatially homogeneous mod- 
els studied previously. For example, in a study of FRW models |l[ it was found that all 



negatively-curved FRW models evolve from the solution corresponding to a global source 
in which the curvature is (negative and) dynamically important (but with a trivial ax- 
ion field) towards the dilaton- moduli- vacuum solutions [^, even in the presence of spatial 
curvature. The physical interpretation of these models, where both the NS-NS two-form 
potential and RR three-form potential are dynamically significant, was discussed in ^1 



with particular emphasis on the fact that the RR field causes the universe to collapse, but 
the NS-NS field has the opposite effect, whereby the interplay between these two fields 
leads to the models undergoing bounces. In the models under investigation here, orbits in 
the full phase space (with ip monotone) approach the dilaton-moduli-vacuum solution on 
the zero- curvature boundary (at late times) and again exhibit a 'bouncing' behaviour; this 
bouncing behaviour is the result of the orbits shadowing orbits in the boundary that are 
constantly being redirected to saddle points of the same or higher stability until it reaches 
a stable equilibrium. 
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